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emissions. NO self absorption is a necessary inclusion of any atmospheric nitric oxide 
analysis stratagem. Profiles of temperature versus altitude and NO column density versus 


altitude for the rocket flight are estimated. 


iil 


TABLE OF CONTENTS 
L. INTRODUCTION oocciiicncs coicccceccceccctnansccccucscsecsessoocsscsusMectaecsscoveuelveeusce: sana ] 
A. THESIS OBJECTIVES. ..........0..0ccscscocecececsedteeucesscssneesss:c sss sttt === eee 3 
B. THESIS QUTLINE. . scce.csessscedeaeetes: 0 ee 3 
Il. BACKGROUND ...............0.cccccccecccescecsceccccccscesesaeseaeesasssegsseceesscueauuueses sss cnsiuniaana— 4 
A. THE ATMOSPHERE. ........ccccccccccecccseccccccccceccccuccocceccsesccececssscee eee sssteiinntainatntniIItInIII=I====——— 4 
d General Description oo... liiccissctsiesosssedtttotc es totttti is ss.s. ttt ee 4 
2. The LOWOSDICTE .........0...c0cc00esadeasseeceeencuieuses s+. sedneeeross0eeneereso<eeggullnnne ie 5 
3. ALMOSPHEVIC ATF QIOW’...........ccccscceccncannanteeseeccccssauseseeseeesapmnsensaesssseyetis «soi 8 
B. MOLECULAR. SPECTRA wx: -csusties ase eee ee ee 8 
Ill. DAYGLOW EMISSION FEATURES ...0....... ccc ccc cc cc cccccccccccceccccccececceceucucescecenees 12 
A. MOLECULAR NITROGEN ...............sccsoccscccsceceeccececareseseccocoeneveityecceesccueccssues «sii 12 
1. Vegdrd-Kaplan Bands ..........iiccccccccccssssssdentusstussesesseseceststiveseessstieseitnss.. tna 12 
2. Lyman-Birge-Hopfield............:isisss.....0iseeeeeees 12 
PSPINTUIR IC OXIDE .......ccccccccecsttettrs occ sc cs etsegiass + « .cacicdien oc <0 elk a cessed se 13 
L. Delta Bands iicccaceiciaviinicckiniwn.Lo ee. eee 13 
2. GONTIIG BOIS ooo. . oc -000c0.ccesvees0eeeesescsuuesse+asscsciekeu cee 13 
3. EpsiloneBand.......... e000” WOR cc ccccccstettrtintte ress oOo nnn etter 14 
IV. THE EXPERIMENT. ...............0....0.0.00.0c00c000UUNDUUINY........0000000. sills. »» nn 16 
A. INTRODUCTION ......c...cccccccceccccocececeucoscsccuccuenceusecuscscccocsccsceseceecsscesscssceltii= III 16 
B. INSTRUMENT DESCRIPTION.............:.cscccocececeseeoscecsceceveecevesoctecscesnccsse cesses se ytiiinnnnnnnnn 16 
C. WAVELENGTH CALIBRATION .........:.02.0ccecescseseesoeccecorcsecvegseeserersseresveses +s see suena iz 
D.DATA COLLECTION ..........0.:cecese++<seacewnncsss+- 200s 0s coun een eer ne nen 18 
V. DATA ANALYSIS ..........0..cceccccccccsscoscccosececcsescssccceecvescosceveseseevesse se seein 19 
A. ALGORITHM TECHNIQUES Sis a wo ES Sic nun nncc onc coc s/QMRNs suas case e'u-cieial ole UOURRIDN Dicieissa 0:0 ua: sicl MIRE re 19 
1. Synthetic Spe Cent, ........ccccccccvseeesovsaesasseeesnuseusssiuniaaiineissadettncca: +: ope nr 19 
2. Fitting Method oocccccccccsoce000s00ssvsaneeasesesse 1) ee re 
B. INTENSITY SCALE FACTOR DETERMINATION ...........ccccccccceccccceccecevcecccccccccceccesceseesecuees 20 
C. NITRIC OXIDE SELF-ABSORPTION ................0.ccecsseccececeosecccecocessccococeracess sess ses usta ZA 
D. NO COLUMN DENSITY DETERMINATION................02.2:0seecceesoeedteceseccacceseeee sess: 21 
E. TEMPERATURE ESTIMATIONS.........cccc0.-..0sescccccoscescosceccecscocectevccsccocece sees: sean Ze 


DUDLEY KNOX LIBRARY 
NAVAL POSTGRADUATE SCHOOL 
MONTEREY CA 93943-51071 


MMMMPIME SUITS .......-0-.-00ccc-200c.0ee0c0c0e0c00000eee-ccssacscceeoscosensteceesecteccessctersnsaceness ee: 23 
A. COMPARISON OF DATA AND SYNTHETIC SPECTRA. .............cccecccccecececccececccsescecsceesceees 23 
Sam CaCO) ONDE ANID. «essen oc cae as scaccsceceace.sdcceccuseeancncccccascscecccecececeeesaccecscceee 24 
MII PBCULAR NITROGEN BSANDS.............0.0ccccceccoscscscccscscsccccescsccosccecscecccscecsescscccsccosseece 25 
mini OXIDE SELF=ABSORPTION ...........-......0-csccerescsscsccsoscoscscocacceccscacceccescecasceccecence 25 
E. NITRIC OXIDE COLUMN DENSITY PROFILE .............ccccccccescecceccecccecceceeceuceccusussuucesensens 26 
PNEIAPERA TURE PROFILE ..........-..-..-...c0-0-ceccaencocsceccosoonrecsconcccescscoseccoccccccssccssccscaccaccsoes 26 
MMEPOPING BLUSIOING............-cccccccescccoccoscscoccscesccsscccescosccsecescescececescecassescesescececsscecescccs De 
RM ey CO NIG S 6... so seis ......0-..-.000c-0----2..ccccccccscscoccscscsccccccsccccscscsccccscscsccccscecs ii 
B. RECOMMENDATIONS FOR FURTHER RESEARCH ...........0cccccccccececcececceccececeucucescuseucececens 24 
VIll, APPENDIX A: FIGURES OF DATA AND SYNTHETIC VERSUS WAVELENGTH ........ 29 
IX. APPENDIX B: NITRIC OXIDE COLUMN DENSITY PROFILE. ...............ccccecceceeceeeees 52 
X. APPENDIX C: TEMPERATURE PROFILE. ............ccccccccccccceucecceccucencccescestecencucescues 54 
XI. APPENDIX D: BEST FIT DATA SCALE FACTORS ...........c.ccccccccccccccececcccececececesees 56 
XU. LIST OF REFERENCES ............... cn ee aa eee 61 
PMMPISU IS IOOGRAPAY ....................0..ccccccccccccsceccccccccscosccscessesceccscccccccecceccscescscscceccccess 62 
XIV. INITIAL DISTRIBUTION LIST ..0.000000000.0oc ccc cccccccccceccccccccccccccccuscccnsceueceescs 63 


ACKNOWLEDGMENTS 
This thesis could not have been completed without exceptional patience of my 
thesis advisor whom I also consider my friend, David D. Cleary, Ph.D. 
My most special thanks to my wife Lict who has endured many evenings of my 
absence to finish this work. I reflect on her as a woman of uncommon beauty and of 


enduring friendship. I will be forever grateful. 


VI 


I. INTRODUCTION 

The electromagnetic properties of the earth's atmosphere are important to many areas of 
science such as weather forecasting, atmospheric instruments, satellite electromagnetic 
surveillance, and communication operations. The upper atmospheric region known as the 
ionosphere contains electron densities which are sufficiently high to form a plasma medium. 
The plasma-like properties of the ionosphere are of specific interest to the military. They 
strongly affect operations of systems such as high frequency (HF) radio communications, over- 
the-horizon radar (OTHR), ballistic missile early warning (e.g., SDI), and the Ground Wave 
Emergency Network (GWEN). For effective prediction of the plasma-like properties of the 
ionosphere an electron density profile is required. 

An electron density profile is a plot of electron density versus altitude for a given time and 
location. In 1986, the Joint Chiefs of Staff priontized the determination of electron density as 
fifth on a list of 43 critical geophysical parameters (The Joint Chiefs of Staff, 1986). The 
electron density varies with time (e.g., day, night, solar cycle), altitude, and geographic 
location. Presently, electron density profiles are determined using approximately 20 ground- 
based ionosonde stations worldwide. This can be very effective locally, but is inherently 
imprecise for a global forecasting model. 

An alternative to ionosonde ground stations is to deduce density profiles of the neutral 
atoms and molecules in the ionosphere from naturally occuring airglow emissions. From these 
neutral density profiles and knowledge of atmospheric photo-chemistry one can infer electron 


and ion density profiles. The inference of electron densities from emission observations is 


performed by relating measurements of natural ultraviolet emissions to the neutral species' 
densities. 

An operational ultraviolet spectrograph that is relatively small and low-powered can easily be 
incorporated into an orbiting tactical communication platform or an optical surveillance satellite 
configuration. It would achieve continuous, passive measurements of the ultraviolet radiation 
emitted by the ionospheric airglow. A low earth orbit (LEQ) based satellite would provide a 
globally complete remote sensing technique of predicting an electron density model of the Earth's 
ionosphere. 

The Naval Research Laboratory (NRL) and The Ionospheric Physics Group at the Naval 
Postgraduate School (NPS) have a joint spectrographic project. The NRL project is a High 
Resolution Airglow and Aurora Spectrograph (HIRAAS) with an observation range of 500 A to 
1500 A. The addition of the HIRAAS/MUSTANG experiment represents a significant progress of 
the global determination of electron density profiles beyond ground-based ionosondes. A sounding 
rocket-launched experiment was first successfully tested in March, 1990 from White Sands Missile 
Range, New Mexico. The data for this thesis is from the second successful launch which was 
conducted in March, 1992. A third rocket experiment will be conducted in March, 1994, from the 
Poker Flat Missile Range in Alaska. The relative high latitude of the Poker Flat Range will enable 
the observation of the effects of auroral electron precipitation on ionospheric photochemistry. 

These sub-orbital rocket experiments will be followed by a remote sensing satellite project. The 
successor to the MUSTANG instrument is the Ionospheric Spectroscopy and Atmosphcric 
Chemistry (ISAAC) instrument which is already constructed. The NPS lonospheric Physics 
Group is prcsently calibrating ISAAC. It will bce launched on the Advanced Rescarch and Global 


Observation Satellite (ARGOS) in late 1995. 


A. Thesis Objectives 

This thesis focuses on the analysis of spectral data obtained from the March, 1992 rocket 
flight. The analysis extends methods used previously on subsets of the data to the observed range, 
1900 A to 3100 A (see e.g., Cleary, 1986; Clayton, 1990; Anderson, 1990; Mack, 1991; Walden, 
1991). All known emissions in this range are used in the fitting routines except OI 2972 A, OII 
2470 A, OIll 2853 A, and NII 2143 A atomic lines. The result of the fitting algorithm is a set of 
relative scale factors as a function of altitude. The relative scale factors are converted into 
absolute terms and plotted against altitude. These absolute plots of intensity and column density 
can be used to infer an electron density profile (see e.g., McCoy, 1985; Meier, 1991). 

Additional goals of this thesis are to confirm the theoretical self absorption mechanisms of the 
NO (0,0), (1,0), (2,0) y band emissions, determine the NO column density profile and estimate a 


temperature profile for the observed altitude range. 


B. Thesis Outline 

This thesis is separated into seven chapters. Chapter II gives general background information 
on the descriptions and designations of the atmosphere and airglow emissions. A basic explanation 
of molecular spectra is also outlined in this section. Chapter III addresses the specific dayglow 
emission characteristics in the wavelength range of interest. Chapter IV presents an overview of 
the experiment. Chapter V gives a full breakdown of the algorithm techniques to generate a 
synthetic model to fit the data. Chapter VI explains the results of the study. Chapter VII 
concludes the document with a summary of findings and suggestions for future study. The 
appendices contain the plots of the final fits to the data as well as an estimated NO column density 


profile and atmospheric temperature profile. 


Il. BACKGROUND 
A. The Atmosphere 
ie Gene Description 
The Earth's atmosphere is designated by distinct, horizontal layers. These layers are 
determined by one of two primary methods. One method designates the layers based upon 
temperature characteristics and the other uses constituent density to classify the atmospheric 
boundaries. Both schemes are shown in Figure II-1 . The figure shows two temperature 


profiles corresponding to low and high solar activity. 
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Figure II-1_ Atmospheric temperature profile. Also shown are atmospheric layers designated 


by temperature compared with corresponding density layers. 


This method of classification designates the troposphere as the region closest to the 
Earth's surface. In the troposphere temperature decreases with altitude, on average from about 


300 K at the surface to 220 K at 10 km. The point at which temperature begins to nse is called 


the tropopause and it determines the start of the next region, the stratosphere. The stratosphere 
is a relatively warm layer due to absorption of solar radiation by ozone (primarily 2000-3000 
A). The top of the stratosphere is marked by the stratopause. The next layer is labeled the 
mesosphere. In the mesosphere, the temperature decreases with altitude until a low of about 
100 K occurs at approximately 85 km. This 1s the mesopause, which marks the beginning of 
the next layer, the thermosphere. Temperature mses steadily in the thermosphere with 
increasing altitude until reaching an equilibrium temperature between 1000 K and 2000 K, 
depending on solar activity. 

Studying the atmosphere in terms of density or chemical composition gives rise to a 
different series of atmospheric layers. The lowest 80 kilometer layer of the atmosphere is 
called the homosphere. Vertical mixing of this region provides a uniform relative density ratio 
between various molecular species. The primary components are molecular nitrogen and 
oxygen and atomic argon in the well known ratio of .78: .21: .009. The remaining 0.001 is 
primarily helium, hydrogen, and carbon dioxide. The layer from 80 km to 100 km is 
designated the barosphere and it is a transition region above which very little vertical mixing 
occurs. Atop the barosphere is the heterosphere from 100 km to approximately 1000 km. In 
this area of the upper atmosphere, the relative densities vary with altitude due to differing 
molecular weights. The density of the atmosphere above the heterosphere is so low that 
individual molecule trajectories are essentially parabolic orbits. This region is labeled the 
exosphere. 

2. The Ionosphere 
The data analysis in this work pertains to the lower heterosphere. This corresponds to 


the ionosphere. Specifically, the first 200 km of the thermosphere where the convective mixing 


is minimal and there exists a positive temperature gradient. The ionosphere is so designated 
due to the ionization of neutral atmospheric constituents into ions and electrons. The 
ionosphere is divided into four layers due to the presence of discrete electron density regions 
and is depicted in Figure II-2 . The layers are labeled D, E, F; and F2. The actual height and 
thickness of each layer varies with solar cycle activity. The D region is the lowest layer and 
typically extends from about 50 km to 90 km. The E region is from 90 km to 130 km. F; 


ordinarily exists from 130 km to 160 km during daytime periods and the F) region lies above 
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Figure II-2 Atmospheric electron density profile. 


The primary mechanism for the development of the ionosphere during the daylight 
hours is photo-ionization. The general chemical equation below for a single, neutral molecule 
M describes this process. 


Mt+hvo > M* +e 


The Av term represents solar photons with sufficient energy for ionization. The ionization rate 
Q,(z) is proportional to: 1) the flux of photons as a function of wavelength /,(z) 2) the 
density of neutral molecules n(z) exposed to the solar flux;; 3) the efficiency with which these 
photons cause ionization. This efficiency is termed the ionization cross-section oA). It is 
wavelength dependent and depends on the ionization potential of the molecule 4. Thus, 
Q(z) = n(z) *1,(2) *o,(A) 

The photon flux decreases in intensity with decreasing altitude and the neutral species 
concentration decreases with increasing altitude. The product of these two quantities produces 
a profile with a discrete layer. This profile is referred to as a Chapman ionization profile. This 
is depicted in Figure II-3. In region (A) of this figure the species concentration is not high 
enough to cause significant photon absorption. Maximum ionization occurs at a distinct 
atmospheric area (B). The ionization rate falls off rapidly at lower altitudes because very few 


ionizing photons are able to penetrate through to the (C) section. 
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Figure II-3 Chapman Profile for 1on/electron production. 


Chapman (1931) was the first to provide a quantitative explanation for the formation 
of an atmospheric ionized layer. This ionization process is effective enough to maintain 
sufficient quantities of ionized species to maintain an ionospheric plasma. The collective 
behavior of the ionospheric plasma can affect electromagnetic waves propagating through the 
ionosphere and is of special interest to the military. 

3. Atmospheric Airglow 

The MUSTANG instrument was designed to measure emissions of the Earth's 
airglow. Airglow is a naturally occurring radiation from the Earth's atmosphere. Dayglow, or 
daytime airglow is primarily driven by photoexcitation. Also contributing to dayglow are 
emissions from photoelectron impacts and photochemical reactions. The intensity of emissions 
resulting from photochemical reactions is negligible during daylight hours compared with the 
intensity from photoexcitation. Photochemical reactions are significant at night since there are 
obviously few direct solar photons. Photochemical reactions are not specifically addressed in 
this analysis because the data was obtained during the daylight period. 

Photoexcitation is the excitation of an atom as a result of the absorption of a photon. 
The atom drops back to a lower energy state by emitting a photon with energy equal to the 


transition. This process will be further described in the next section on molecular emissions. 


B. Molecular Spectra 

Molecular emissions may occur as the result of electronic state changes, but they may also 
occur from purely vibrational and rotational changes of state. Figure I-4 illustrates an 
electronic transition. Each of these types of transitions has typical energy gaps, with electronic 


being on the order of 10 eV, vibrational of .1 eV, and rotational of .001 eV. 


As seen in the figure, for any given electronic energy level there are many vibrational 


energy levels labeled v’ for the excited electronic state or v" for the lower or ground electronic 


state. Similarly, for a given change in vibrational level there are many possible initial and final 


rotational levels. As a result, molecular emissions from one vibrational to another are 


composed of many closely spaced lines. These emissions are called molecular bands. 


The principal photon-molecule interactions are shown in Figure H-5. For elastic 


scattering, an energy transfer between the photon and the molecule does not take place, even 


though the photon may change direction. With inelastic scattering, the incident photon not only 


changes direction, but imparts energy to the molecule. In this case, the molecule is excited to a 
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Figure Il-4 Molecular energy versus internuclear distance. 


higher energy state. If the energy difference between states 1s equal to a level AE, then the 
deflected photon will propagate with an energy AE less than the initial photon’s incident 
energy. Another photon with this energy AE will be emitted when the molecule retums to the 
ground state. Resonance scattering 1s initiated when the incident photon is absorbed, raising 
the molecule to an excited state. This is followed by a re-emission when a photon with 
wavelength equal to the incident photon is emitted. Fluorescence scattering occurs when the 
excited molecule undergoes several downward tcnsitions in returning to the ground state. 
Each transition generates a separate photon. An excited molecule can also collide with another 
particle, thus returning to the ground state without emitting a photon. This 1s called collisional 
de-excitation or quenching. 

Similar to photon-molecule collisions, electron-molecule collisions may be elastic or 
inelastic. With inelastic collisions, the molecule is excited to a higher energy level and the 
electron loses kinetic energy. Photon emissions occur when the molecule drops from the 
excited state to the ground state. The primary interactions considered for this study are photon- 


molecule inelastic scattering and flourescense scattering. 
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Figure II-5 Molecular processes induced by incident radiation. 
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Ill. DAYGLOW EMISSION FEATURES 

A. Molecular Nitrogen 

Figure III-1 is an N, energy level diagram including the two transitional bands of interest 
in the 1800-3400 A range. The states are indicated by molecular notation, and the transitions 
are identified by name. Since these are forbidden transitions, the excitation is from electron- 
molecule inelastic scattering. 

1. Vegard-Kaplan Bands 

The Vegard-Kaplan system 1s so designated because Vegard discovered the system 

in 1932, and Kaplan determined the specific wavelengths involved a few years later. The 


emissions range from 1250 A to 5325 A and represent the Ab2,*>X!2,* transition. In this 


transition, the total transformation is +—>+, the electronic inversion is u>g, AZ= -1, and AA=0. 
From the molecule selection rules, this is a forbidden transition, yet it has been observed as an 
electric dipole emission. A thorough discussion of the band system is given in Meier (1990). 

A-state nitrogen molecules are produced by cascading downward from the B and C 
states or as a result of photoelectron collisions. Because the relationship of emission bands to 
photoelectron flux is complicated by cascading N2 molecules and quenching effects, the V-K 
band is not generally used to infer photoelectron flux. 


2. Lyman-Birge-Hopfield 
The Lyman-Birge-Hopfield system is the name given to the ATI, >X1Zig* transition 
of N,. The emission bands arising from this transition extend from 1160 A to 3020 A and are 


due to both electric quadrupole and magnetic dipole transitions. Unlike the V-K system, the 


upper state of the LBH system can only be populated by direct photoelectron collision 
excitation. This system can, therefore, be directly used to infer the photoelectron flux. The 
LBH bands found in the MUSTANG data are observed in the wavelength region from 1800 A 
to 2400 A. An intensive investigation of the LBH emissions observed during the 1990 flight 


was performed by Mack (1991). 


B. Nitric Oxide 
Nitric oxide (NO) emissions dominate much of the dayglow data. There are four primary 
band systems in the range 1800 A-3400 A: 7, 5, B, and. These NO transitions are generated 
by photoexcitation. The B bands are several orders of magnitude less intense than the other 
emissions, so they have been neglected in this analysis. Figure II-2 is an NO energy level 
diagram including the three transitional bands of interest in the 1800 A-3400 A range. 
1. Delta Bands 
The C2II—>X2I transition band system, designated 5, ranges from 1600 A to 4300 A 
(Barth, 1965). Using standard notation, the upper electronic state's vibrational level is 
designated v', and the lower state vibrational level, v". The 5 transition includes vibrational 
levels corresponding to v'=0,1,2,3,4 and v"=0,1,...23. Nitric oxide will dissociate at v' > 0 for 
CTI, so v' is limited to v= 0 only. 
2. Gamma Bands 
The y bands of nitric oxide extend from 1850 A to 6150 A. The emissions arise from 
the A?2’—>X2I transition, and the ground state has v"= 0,...23 as previously stated. The 
dissociation energy for nitric oxide is lower than the energy of the v'= 4 level in the A state, so 


only transitions arising from v= 0,1,2,3 contribute to the synthetic y spectrum. 
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3. Epsilon Bands 
The nitric oxide € system is the result of molecules transitioning from D?= err 
The wavelengths of emissions range from 1800 A to 2400 A. D-state molecules populate 


vibrational levels v'= 0,1,2. 
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IV. THE EXPERIMENT 
A. Introduction 
The middle ultraviolet spectra analyzed in this thesis were obtained from the NPS 
MUSTANG spectrograph flown on a NASA sounding rocket launched on March 19, 1992, 
from White Sands Missile Range, New Mexico. The MUSTANG instrument was launched on 
a Terrier boosted Black Brant vehicle, which carried the payload to an apogee altitude of 
320km. The spectrograph field-of-view was pointed in the anti-solar direction with an 


observation-zenith-angle of 100° on the up-leg and 90° on the down-leg. 


B. Instrument Description 
The MUSTANG (Middle Ultraviolet Spectrograph) shown in Figure IV-1, covers a 


wavelength range from 1850A to 3400A. It is a modified 1/8 m Ebert-Fastie spectrograph 


1/8th m Telescope Mirror 


Spectrograph Entrance Slit 





] 1/8th m Ebert Mirror 


ITT Image Intensifier hy 


Hamamatsu Image Sensor 


Figure IV-1 Schematic diagram of the MUSTANG instrument. 
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with a 1/8 m telescope constructed by Research Support Instruments, Inc. The entrance slit is 
5 mm by 140 ym which yields a field-of-view of 2.3° by 0.07°. A 10 A wavelength resolution 
is provided by a HYPERFINE (Inc.) 1200 lines/mm grating ruling. The detector assembly | 
consists of an ultraviolet-to-visible ITT image intensifier and a Hamamatsu linear array 
photodiode detector. The image intensifier uses a quartz input window, a CsTe photocathode, 
and a fiber optic output window. The detector is an array of 512 photodiodes of area 5 mm by 
50 um giving analog voltage outputs proportional to intensity using a fixed 50 msec integration 
period. This voltage output is fed serially to an A/D converter and then to rocket telemetry via 
a FIFO shift register. 

Pre-launch calibration/sensitivity was obtained using standard techniques with a 
Deuterium lamp for the wavelength range 1800 A-2800 A and a FEL tungsten lamp for the 
range 2400 A-3400 A. The calibration curves obtained from the two lamps agree to within 
10%. The MUSTANG instrument experienced no significant degradation in sensitivity (less 


than 5%) during the launch sequence based on post-flight calibration. 


C. Wavelength Calibration 

The wavelength calibration was determined using the emission spectra of Mercury and 
Platinum. The platinum lamp provides 25 spectral emissions in the wavelength region of 
interest. The design of the MUSTANG instrument is optimized for a bandwidth from 1850 A 
to 3400 A and a dispersion of 3.13 A/pixel. 

The pixel positions on which the selected platinum emissions fell were determined by a 
rough approximation ( + 0.5 A) of the centroid of each spectral line. The wavelengths of the 


25 emissions were plotted versus pixel position. Previous works (Anderson, 1990: 
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Chase,1992) used a linear best fit calibration procedure. It was decided to attempt a 


polynomial best fit for this analysis. The polynomial fit was determined to be: 
A =[1848.270+ 3.188N — (1374X107)N?| A; 


where JN ts the pixel number. 


D. Data Collection 

Over 8000 spectra were collected by the MUSTANG instrument during the rocket 
experiment. These spectra were averaged into 65 altitude bins, disregarding a few spectra 
which were contaminated by noise spikes from the instrument's electronics. Each of the 65 
averaged spectra covers a 5 km layer that is labeled by its central altitude. The adjusted 
spectra were then divided by instrument sensitivity. Each spectrum contained 512 intensity 
values. The first two values were set to zero, because the uncertainty in the instrument 


sensitivity at the lower wavelength endpoint is extremely large. 
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Vv. DATA ANALYSIS 
A. Algorithm Techniques 
1. Synthetic Spectra 

Computer synthetic spectra written by Cleary (1986) were used for the generation of 
Nitric Oxide 5, ©, and y and the Nitrogen Lyman-Birge-Hopfield band emissions. The 
Nitrogen Vegard-Kaplan synthetic model was written by Siskind and Barth (1987). All of the 
computer code used in this study were written in Interactive Data Language (IDL*), Version 
3.1, designed by Research Systems, Inc. The inputs for the synthetic spectrum algorithms 
were temperature and the calibrated wavelength range of the MUSTANG instrument. The 
output is an array of 512 intensities corresponding to the instruments wavelength range. 
Before comparing the synthetic spectra with the data, the synthetic spectra were convolved 
with MUSTANG’s instrument response function. The MUSTANG instrument response was 
estimated to be an 11-element symmetric slit function (Clayton, 1990; Anderson, 1990; Mack, 
1991). 

2. Fitting Method 

The technique used to fit the synthetic spectra to the data was an iterative process of 
modifying a linear combination of each band emission system and adjusting separate synthetic 
and quantum parameters of each emission as necessary. These adjustable parameters are the: 
1) Intensity Scale Factor; 2) Molecular Oscillator Strength; and 3) the Franck-Condon Factor. 

The intensity scale factor is a relative weight assigned to each band system. It is a 


best “initial guess” taking into account the combined contributions of other surrounding 
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emissions. For example, the entire N2 LBH model was assigned a specific scale factor and 
then summed with all other contributing spectra. 

The molecular oscillator strength 1s a quantum mechanical component for the 
synthetic models which quantifies the probability of a molecule absorbing a photon in the 
electronic ground state and transitioning to a specific vibrational level (v') in the excited 
electronic state. This mechanism of molecular photon absorption is discussed in the 
background section. For instance, in the N2 VK synthetic band system, at a temperature of 
1125 K, all v'= 3 bands were assigned a scale factor of .7 in order to decrease the intensity by 
30% for a best fit. 

The Franck-Condon factor is a quantum mechanical element which quantifies the 
probability of a synthetic emission occurring from a specific vibrational level (v') in the excited 
electronic state to a specific vibrational level (v") in a lower or ground electronic state. For 
example, within the NO y synthetic band a the y (2,3) emission was observed to be 
underestimated by 50% but the y (2,6) was overestimated by 35% and both were assigned a 
specific scale factor to correct the discrepency. 

The adjustable synthetic and quantum parameters designated as the intensity scale 
factor, molecular oscillator strength, and the Franck-Condon factor will be referred to as the 


“SQ parameters” for the rest of this work. 


B. Intensity Scale Factor Determination 

Intensity scale factors were determined by visual inspection for “goodness of fit”. An 
analysis strategy was deveioped to take advantage of various attributes of NO and N>2 
molecular spectra. For example, the N2 LBH and NO e band systems are more prominent 


from 1900 A to 2200 A while the N, VK and NO y band systems are more dominant in the 
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2200 A to 3100 A range. A time intensive iterative approach was conducted to identify each 
synthetic band system which had a discrepency with the data. An initial attempt was made to 
assign one scale factor per band system. However, it was necessary to modify individual 
emissions. The SQ parameters of these emissions were modified to achieve a best fit synthetic 


spectrum to the rocket data. 


C. Nitric Oxide Self-Absorption 

A self absorption model using Holstein transmission functions, developed by Eparvier and 
Barth (1992), was applied to the NO (0,0), (1,0), (2,0) y resonance band emissions. The 
generated synthetic models were then added to other contributing spectra. If the intensities had 


a discrepency to the data, their respective SQ parameters were modified for a best fit. 


D. NO Column Density Determination 
The NO column density profile determination was performed in tandem with the 
temperature profile and NO intensity scale factor estimations. The NO column density (pci) 


for each altitude bin was calculated from the following relationship: 








atoms 
Pos = nr (S28 
we -sn{A| = NO intensi le factor (* 
heal pixel = intensity scale factor ii 
R=1X 10° (photons) pref annie ce) 
cm’ * sec photon 


The « term is the MUSTANG instrument pixel spacing. The dependence on instrument 


sensitivity has been removed from the rocket data spectra to yield absolute spectral intensity 
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R 
(n) in units of Rayleighs per Angstrom (} A Rayleigh (R) is defined 


; { photons . | | 
as 1X10 eee And the ¢ term is a conversion factor representing one photon 


emitted from each NO atom and detected by an individual pixel of the MUSTANG instrument 


per unit time. 


E. Temperature Estimations 

The temperature of each altitude bin was estimated by isolating an observed molecular 
band relatively free of other contaminating spectra. The standard emission used below 160 km 
was the NO y (0,1) at 2365 A. Above 160 km the Nz VK system begins to contaminate the 
NO y (0,1). The N2 VK (0,5) at 2605 A was used for estimating a temperature profile for all 


of the upleg and downleg altitude bins above 160 km. 
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VI. RESULTS 
A. Comparison of Data and Synthetic Spectra 

The averaged spectra examined here are for the altitude data bins 110, 160, and 195 km. 
The primary emission features observed in this MUV spectral region of 1900 A to 3100 A 
include the Nj Lyman-Birge-Hopfield (LBH), Nz Vegard-Kaplan (VK), and the NO Gamma 
(y), Delta (6), and Epsilon (€) bands. All known emissions in this range are used in the fitting 
routines except OI 2972 A, Oll 2470 A, OMI 2853 A, and NII 2143 A atomic lines. 

A comparison of the observed spectra with the synthetic spectrum model in the spectrum 
for all three of the altitude ranges is shown in Figures A-1 to A-22 in Appendix A. In order to 
illustrate the basic technique and results of the fitting technique employed for this work, the 
160 km data set was used as a sample (Figures A-1 to A-10). The 160 km data set was chosen 


because the intensities of the NO and N> band emission systems were approximately equal. 


Figure A-1 illustrates the 160 km NO « (v'=0) emissions for the range 1900 A to 2130 A. 
The N2 LBH (v'=3) emissions in the same wavelength range are shown in Figure A-2. Both of 
these synthetic spectra are given to show the individual contributing emissions of the five band 
systems. To show a complete band system, the sum of all the N2 LBH emissions is shown in 
Figure A-3. Figure A-4 presents the summation of the N2 LBH and NO e band systems. 
Continuing with the step by step additions of the other three band systems, the result is shown 
in Figure A-5. 

The individual emissions were used to identified discrepancies between the data and 


synthetic spectra. This was accomplished by overplotting several emissions in the wavelength 
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range in question. Once an element was identified as the erroneous contributor (e.g., exceeded 
or underestimated observed intensities), its respective SQ parameters were modified then 
resummed for a best fit. 

Figure A-5 reveals several interesting features of the 160km data set. Note that only 
minimum intensity scale factors are applied. One can visually identify the synthetic spectrum 
exceeding the observed data at approximately 1910 A, 1980 A, 2050 A, and 2095 A. The 
synthetic spectrum can also be seen lagging at 2030 A and 2025 A. 

Figure A-6 presents the data versus wavelength compared with the theoretical, but 
modified oscillator strengths to v' emissions have also been applied. Notice that the model 
underestimated the observed intensity at 1910 A and 1930 A. To correct this, Figure A-7 
shows how modifying the Frank-Condon factors for each erroneous (v', v") emission adjusted 
the theoretical fit closer to the real data. Figures A-8, A-9, and A-10 repeats this process and 
displays the entire 160 km data set for 1900 A to 3100 A. In addition, all of the altitude data 
bins at 110, 160, and 195 km are illustrated in Figures A-11 to A-22. 

The resulting quantum molecular oscillator strengths (v') and Franck-Condon factors (v") 
to achieve a best fit for each band emission for all three altitude data sets are listed in Tables 1 


through 3 in the Appendix section. 


B. Nitric Oxide Bands 

The nitric oxide synthetic spectra were substantially contaminated by molecular nitrogen at 
160 km. The y (0,1) is an intensity standard at lower altitudes but is interfered with by 
surrounding N2 Vegard Kaplan emissions. 

A successful attempt was made to standardize the NO Gamma (y), Delta (5), and Epsilon 


(c) bands by weighing all three systems by the same intensity scale factor. This was 
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accomplished by modifying key values of the SQ parameters for each band emission. The NO 


Delta (5), and Epsilon (¢) bands were the most difficult to model due to contaminating N2 


Lyman-Birge-Hopfield (LBH) emissions in the 1900 A to 2120 A range. The NO Gamma (y) 


band system was also contaminated with Ny VK emissions in the 2200 A to 3100 A range, but 


to a lesser degree. 


C. Molecular Nitrogen Bands 

At 160 km the Vegard Kaplan bands begin to dominate over all other emissions with the 
exception of the y (1,0). Initial analysis indicated two discrepencies of the synthetic spectra 
concerning these bands. First, a negative 3.13 A shift (approximately 1 pixel resolution) of the 
VK synthetic versus data spectra. To correct this discrepency, a modification to the VK 
synthetic algorithm (VKTERM.PRO) of the initial molecular constants was performed. The 
TO molecular constant was change from 50203.6 to 50169.0 and the WXEE molecular 
constant was changed from 14.32 to 11.45. And the second discrepency was the observed 
relative intensities of the VK (0,5), (0,6), (0,7) bands. The synthetic predicted a constant ratio 
of these three emissions over all temperature inputs. However, the observed data showed 
varying scale factors which were matched by modifying the synthetic model and have been 


recorded in Appendix D, Tables B-1 through B-3. 


D. Nitric Oxide Self-Absorption 

A self absorption model was made of the theoretical self absorption versus non-self 
absorption mechanisms of the NO (0,0), (1,0), (2,0) y resonance band emissions. The Holstein 
transmission functions applied for this model generated mixed results. The NO (0;0) y 


resonance band was consistently overestimated at all three altitudes analyzed. The NO (1,0), 


25 


(2,0), y resonance bands also had to have modifications to their respective SQ parameters, but 
the results did not indicate any consistency. Figure A-11 through A-14 includes the NO y 


resonance bands at 110 km. 


E. Nitric Oxide Column Density Profile 
The nitric oxide density profile inferred from the data is determined from all of the 
recorded altitude bins and is shown in Appendix B, Figure B-1 and Table B-1. These 


estimated densities are consistent with other observations during simular solar cycle intensities. 


F. Temperature Profile 
A temperature profile estimated from the data is determined from all of the recorded 
altitude bins and is shown in Appendix C, Figure C-1 and Table C-1. These estimated 


temperatures are consistent with other observations during similar solar cycle intensities. 
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VII. CONCLUSIONS 
A. Summary of Findings 
In this work, spectra from the Earth’s ionospheric dayglow obtained from a sounding 
rocket-launched experiment were analyzed. Computer generated spectra were created for the 


NO y, 5, and ¢ bands and the Nj LBH and VK bands. These synthetic spectra were fitted to 


the observed dayglow data. An overall synthetic fit agreed well with the observed spectra. 
Intensity scaling factors, molecular oscillator strengths and Franck-Condon factors were used 
and modified to obtain a best fit. Synthetic emission bands that were observed to be in 
discrepency with the observed data were visually identified and corrected to obtain a best 
possible fit. 

The NO intensity scaling factors were used to infer a NO column density profile. A 
temperature profile was estimated from the synthetic spectra which best fit the data. The 
synthetic fit also confirmed the theoretical self absorption versus non-self absorption processes 
of the NO (0,0), (1,0), (2,0) y resonance band emissions. NO self absorption is a necessary 


consideration of any atmospheric nitric oxide analysis scheme. 


B. Recommendations For Further Research 

The NPS MUSTANG will be part of an operational remote sensing spacecraft. It will be 
launched on the Advanced Research and Global Observation Satellite (ARGOS) in late 1995. 
Considering the time intensive nature of a visual analysis by comparing the observed data with 
synthetic spectra, a high speed central processor is strongly recommended. Once an intensive 


catalog of synthetic emissions (which are in discrepency with observed airglow) is complete, 


Zl 


appropriate modifications to the fitting programs would add to the reliability of MUSTANG as 


a surveyor of ionosperic parameters for national defense agencies. 
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VIII. APPENDIX A 


Appendix A contains the figures of rocket data versus wavelength and comparing the 
data with the generated theoretical spectrums as referred to in the results section. 
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IX. APPENDIX B 
Appendix B contains Figure B-1 which shows the NO column density profile versus 


altitude and Table B-1 which lists each data point of the NO column density profile as referred to in 
the results section. 
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xX. APPENDIX C 


Appendix C contains Figure C-1 which shows the temperature profile versus altitude and 
Table C-1 which lists each data point of the temperature profile as referred to in the results section. 
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XI. APPENDIX D 


Appendix D contains the listings for minimum scale factors, modified oscillator strengths, 
and modified Franck-Condon factors (Best Fits) corresponding to the data bins 110, 160, and 
195km. and are shown in Tables D-1 thru D-4. 

Table D-1 lists the intensity scale factors for the Nz LBH and VK synthetic spectra, and 


the NO y, 5, and € synthetic spectra. 

Tables D-2, D-3, and D-4 contain the fractions necessary to modify each band emission for 
a best fit. Ifa fraction is entered for a v' row, it is equivalent to applying that percentage to the 
molecular oscillator strength. If a fraction is entered for a v" row, it is equivalent to applying that 
percentage to the molecular oscillator strength and the Franck-Condon factor. If any entry is 
blank, assume the factor 1s 1.0. 
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Table D-1 





Synthetic Scale Factors 






1125 

7.59E+03 1.46E+04 1.59E+04 
4.46E+04 2.50E+04 1.46E+04 
3.51E+09 5.46E+08 2.79E+08 


_NO Column Density: 9.08E+15 zz ES 5.13E+14 


Synthetic Self Absorption Scale Factors 
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Table D-2 Altitude: 110 km Temp: 200 K 

vy _L | VK LBH]| GAM] EPS| DEL| |v] v | VK] LBA) CAM] BPS] DEL 
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| o4 | tT le 
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Pe foe |e eae eee 


58 


Table D-3 Altitude: 160 km Temp: 650 K 
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